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oo Abstract 
O 

We derive a class of formulae relating moments of B — > X u tv to B — > X s ^j in the shape function 
region, where m 2 x ~ i^Aqcd- We also derive an analogous class of formulae involving the decay 
B — > A s £ + ^~. These results incorporate power corrections, but are independent of 

Q leading and subleading hadronic shape functions. Consequently, they enable one to determine 

\Q \V u b\/\VtbVf S \ to subleading order in a model- independent way. 
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I. INTRODUCTION 



The study of decays of the B meson allows us to probe QCD and flavour physics. The 
program's goals include, on the one hand, precision measurements of Standard Model pa- 
rameters and, on the other hand, searches for new physics. Short-distance physics is encoded 
in Wilson coefficients of local operators. By comparing measurements of these coefficients 
with theoretical predictions, signals of new physics may be found. High sensitivity to new 
physics is provided by the so-called rare decays, namely those channels involving flavour- 
changing neutral currents, since they do not occur at tree level in the Standard Model. 
Measurements of the inclusive rare process B — > X s ^y [1, 2, 3, 4, 5] have provided significant 
constraints on extensions to the Standard Model. The more complicated decay B — > X s £ + £~ 
is complementary to B — > X s 7, as its effective Hamiltonian includes two extra operators. 
Moreover, additional observables are available, such as the q 2 spectrum and the forward- 
backward asymmetry, which have been the focus of much work. Recently, it was noted that 
an angular decomposition provides a third observable sensitive to a different combination 
of Wilson coefficients [6]. Belle and BABAR have already made initial measurements of 
B -> X s £ + £~ [7, 8]. 

Precision measurements also provide determinations of elements of the CKM matrix or, 
equivalently, the angles and sides of the unitarity triangle. By overconstraining these, the 
flavour structure of the Standard Model is subjected to rigorous examination. For the decay 
B — > X c £u, experimental and theoretical uncertainties are under control, and consequently 
\V cb \ is one of the best-determined elements of the CKM matrix. From B — > X u £u, we can 
also determine \V ub \ [9, 10, 11, 12]. 

However, inclusive B decays often require a trade-off between theoretical and experi- 
mental difficulty: if phase-space cuts are necessary experimentally, then the spectra will be 
less inclusive and the corresponding theory more complicated. In this respect, B — > X c £v 
and B — > X u £v are markedly different. The former is sufficiently inclusive to enable the 
use of a local operator product expansion (OPE) [13], in which non-perturbative correc- 
tions appear as an expansion in inverse powers of m h . This formalism has been calculated 
to order l/mf [14] (and recently to order l/mf [15]), with the relevant non-perturbative 
matrix elements defined via the Heavy Quark Effective Theory (HQET) [16, 17, 18]. In 
contrast, in B — > X u £v experimental cuts (e.g. cuts on Ei or m 2 x ) are required in order to 
eliminate the dominant b — > c background. In many cases, we are restricted to a region in 
which m 2 x ~ to&Aqcd and the local OPE breaks down. In this so-called endpoint or shape 
function region [19], the set of outgoing hadronic states becomes jet-like and the relevant 
degrees of freedom are collinear and ultrasoft modes. The Soft-Collinear Effective Theory 
(SCET) [20, 21, 22, 23] is then a powerful theoretical method. 

Similarly, B — > X s 7 measurements employ a cut on the photon energy. In Refs. [24, 25] 
it was shown that the shape function region is also relevant for B — > X s £ + £~ . Here, cuts 
are made in the dileptonic mass spectrum to remove the largest cc resonances, namely the 
J/^f and These leave two perturbative windows, the low-g 2 and high-g 2 regions. At low 
q 2 , where the rate is higher, an additional cut is needed: a hadronic invariant-mass cut is 
imposed in order to eliminate the background b — > c(— > s£ + v)£^v. 

At leading order (LO) in Aq CD /mfe, decay rates now depend upon a non-perturbative, 
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and hence analytically incalculable, shape function. However, this function is process- 
independent and appears in both B — > X u £v and B — > X s 7, for example. One can thus 
measure the leading-order shape function from the photon energy spectrum of B — ■» X s 7 
and use the result in the B — > X u £z/ spectrum, or, more directly, express the semileptonic 
rate in terms of the radiative rate instead of the shape function [26, 27, 28, 29]. In this way, 
model dependence can be avoided in the determination of \V u b\. 

At subleading order, the situation is far more complicated, with several universal shape 
functions occurring in different combinations [30, 31, 32, 33, 34, 35]. In this paper, we 
construct combinations of shape-function-dependent decay rates that are protected from 
non-perturbative effects to second order in the power expansion. Through this procedure, 
we obtain formulae for | V u b\/\ VtbV t *\ that are free from the hadronic uncertainties arising 
from the leading and subleading shape functions. This method uses moments of the fully 
differential decay spectra of B — > X u £v and B — > X s 7 (and, optionally, B — > X s £ + £~). 

The rest of this paper is organized as follows. In Sec. 2, together with Appendices A 
and B, we present the basic formalism needed for our work. This includes power corrections 
for the triply differential decay spectra of the semileptonic processes and the photon energy 
spectrum of B — > X s 7. In Sec. 3, we derive and discuss our results, eliminating shape 
functions from expressions for \V u b\ at next-to- leading order (NLO). We conclude in Sec. 4. 



II. FORMALISM 



In this section, we briefly review the formalism and results from Refs. [24, 32, 36] that 
we shall use in this paper (see these references for further details). 

The inclusive decay rate for B — > X u £v (B — > X s 7) is proportional to W^L^, where 
is the leptonic (photonic) tensor and is the hadronic tensor, which can be written 

as 

W v = ^ B Y,^ 5 ^-q-px){B\Jl\X){X\J v \B) 

= -Q^Wi + VpVvWz + ie^ aP v a q^W 3 + q^q v W A + (v^q v + v u q^W 5 . (1) 

Here, is the velocity of the B meson and q^ is the £v (7) momentum. We use the 
hadronic current J (e.g. J™ = u^^Pijo for B — > X u £v) and relativistic normalization for 
the \B) states. Similarly, the inclusive decay rate for B — > X s £ + (p + )£~(p~) is proportional 
to (W^Lf + W*L%), where L% R) = 2 [p^p v _ + ptp u + - g"" p+-p- T 'p+aP-p] and 

can be defined analogously to Eq. (1), in terms of a current J L ( R ) [37]. 
Contracting with W^ v and neglecting the mass of the leptons give the differential 

decay rates 

8E, 
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q * W u + (2E^E + - (fl2)W? + q\E. - E+)W™ 6(4E_E + - q 2 ) , 
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for B -> X s7 , 5 -> XJz/ and B -> respectively, where W^ 2) = W^ 2) + Wj 2) , 

IVjf = W% — Wi = Wi(q 2 ,v-q) and the normalization factors are 

G 2 F m% 2 M2|r» eff (°)/' r r, r» _ G 2 F m% 2 



32tt 4 ''""«" umi uv *" j ' ' v w/l ' ' " I92tt :{ 
192tt 3 16tt 2 |nbV/ ' sl ' 



r « ^r F m B u. em 2 



In SCET, it is natural to use light-cone coordinates, where we introduce vectors n and n 
such that n 2 = n 2 = and n-n = 2. A four-vector then has components (p + ,p~ ,p±) = 
(n-p,n-p,pj_). In the region of interest, the set of hadronic states X is jet-like, i.e. p\ <C p%- 
For convenience we define the dimensionless variables 

2E 1 2E t _ n-px n-p x ... 

%h = — , x h = — , y H = , u h = • (4) 

ttib rriB ttib rris 

In terms of these variables, the decay rates are 

^ =K ^{wl-Wl- mBX lWl}, (5) 



1 d 3 T u 
F% dx H dy H du H 

d 3 T u 



24m B (y H -u H ){(l-u H )(l-y H )W^+ 1 -(l-x H -u H )(x H +y H ~l)W^ 

+ n ^(l-u H )(l-y H )(2x H +u H +y H -2)WZ} , 

i d 6 T w ( 1 

Tf dx H dy H du H = 2AmB ^H-UH){(l-UH)(l-y H )W^ + -(l-x H -u H )(x H +y H -l)W? 

l-u H )(l-y H )(2x H +u H +y H -2)W* e } , 
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where Wi = Wi(uH,yn)- The full phase-space limits are given in Table II of Ref. [32]. 

The optical theorem relates the Wi to forward-scattering amplitudes, which can be cal- 
culated by taking time-ordered products of currents. An important part of the analysis is 
the separation of short- and long-distance contributions. The results, known as factorization 
theorems, may be written schematically in the form 

dT = H x J ® / , 

where <g> denotes a convolution. The hard (H) and jet (J) functions encode perturbative 
corrections that appear at two different scales, /if, ~ m b and fa ~ ^/rn b Aq CD respectively, 
whereas the shape function (/) represents non-perturbative physics. 

SCET involves a power expansion in the small parameter A = a/Aqcd/^6- At leading 
order in A, rates depend on one shape function, which we denote by 

fPx 

W^ 0) = h l (p- x ,m b , f i) / dk + jW(pj c k + ,rifW(k + + A-p+,ri, (6) 

Jo 

where A = m# — rrib + (Ai + 3A2) / (2rrib) + . . .. The first subleading shape functions occur at 
order A 2 and we denote these by /q%, f^l and f^. These are common to the three decays, 
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but appear in different combinations, and are convoluted with jet functions ^J^°\ J^" 2 ^ and 
J^^\ respectively, as shown in Eq. (B8). Note that we also have uh/vh ~ A 2 in the shape 
function region. 

The shape functions are given by 5-meson matrix elements of non-local ultrasoft oper- 
ators. The definitions used here follow Ref. [32] and are included in Appendix A. At tree 
level, the jet functions are 

5(kt) - 5(kt) 



(7) 



j(~ 4 '(^)=4 TOs ( ft ) 
At one-loop order, we have 



s(kt) 6(k*) m) 



+ 



+ 



(ktm) (Km) (Km) 



- 7r 2 5(kt)5(k+)5(kZ) 



1 + ^M^ f 21n2 ^_ 3 ln^ + 7- 7 r 2 ) 



(8) 



^ a s {jj)C F 



4:71 



where uj = n-p is the large partonic momentum. 
For convenience we define 



F(p + ,p~ 



dk + J(°\p-,k + ,n) f°\k+ + A-p + ,ti) 



(9) 



o 



/f(A-P + ) 



2m b 

(2) 



Ai + 3 A 2 
2m i 



/ (0)/ (A-p + ), 



FiAp + ) = /iJ'(A- P + ), 

where a prime denotes a derivative, as well as 



F,M)= dkfdkt 



5(kf) - S(kt) 



F 5i6 (p+) = Jdkfdk+ dk+ 



5{kt) 



/& ) (fc;+A- P +), 

8{kt) 5{kt) 



(10) 



(ktm) (Km) (Km) 



-n 2 5(kt)5(kt)5(kl 



xf$$(k+ + A-p + ). 

If we use the tree-level expression for J(°\ then F(p + ,p~) = F(p + ) is a function of p + only. 
Then, for B — > X s 7, the rate dT s /dx] I in the endpoint region is [32] 1 



1 dT s 



V s dr 1 



Qj TT T 



m B (C^) 2 [1 - 3(1 -xl)] F(m B (l-xl), m B ) 

+ [m B (l-x] I )-A]F(m B (l-x' r H )) + F 2 (m B (l-x] { )) 
-F 3 (m B (l-x' r H ))+F A (m B (l-x' r H ))-8na 8 ( t i i ) F°(m B (l-xD) , 



(11) 



1 This includes Or - Or and Or - 02 contributions only. In Ref. [38] subleading corrections from Or - 
08 are studied and estimated to contribute between —0.3% and —3% to the total flavour-averaged decay 
rate. We do not consider such corrections in this work. 
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where 1 — x c H ~ A 2 and 

C («) = 1 + A l(mt ,,)-^^g + 6}, 

^ (m " s) = C° l " »(m t ) {S C7 ' ff " )( ""' ) +E C ?' <0) (">»)'-*(e)} ■ < 12 > 

The triply differential decay rate for B — > X u £z/ at NLO [32] is obtained by substituting the 
listed in Appendix B into Eq. (5). At tree level, this becomes 

1 d 3 T u f 

™ 3 3^—} — = 6(1 - u H ){x H + Uh - l)s 2m B (2 - x H - Vh - u H )F{m B u H ) 

1 ft dy H du H I 

- (y 2 H ~ (2-x H )y H + 2(1- x H ) - u H {2-x H -u H )) Fi(m B u H ) 

Vh-u h 

2 

+ —ri Avh ~ ( 2 - x H)y 2 H - {^-u H ){xH+u H )yH 

yH{yH-UH) 

+2(x H +2y H +u H -l)) F 2 {m B u H ) 

2 

+ —(x H + y H + u H -2)F 3 (m B u H ) 
Vh 

2 

- — — r (y 2 H -(2-x H )y H +2(l-x H )-u H (2-x H -u H )) F 4 (m B u H ) 

yH{yH-u H > 

4 

- —rz r(l - yH){x H + VH-1) ^na s {^i)F^(m B u H ) 

yniyH-UH) 

4 1 

+ —- r(l - u H ){l - x H - u H ) ^TTa s (iii)FQ(m B u H ) > ■ (13) 

yH{yH-u H ) ) 

Note that we can use the relation [30] 

= 2{K-m B u H )F{m B u H ) + 0{\ A ) (14) 

to eliminate Fi(m B %), as was done in Eq. (11). 

The triply differential decay rate for B — > X s l + l~ was calculated in Refs. [24, 36]. The 
appearing in Eq. (5) are also listed in Appendix B. 



III. \V ub \ AT NLO 

A. Relations between B — > X u tl> and B — > X s 7 

Consider first the process 5 — > X u £z/. We wish to isolate or eliminate the subleading 
shape functions that appear in the rates. In the following, we shall work at tree level. 
Inspection of Eqs. (B2) and (B9) shows that the shape functions appear in the hadronic 
structure functions W\ to W 3 in only two combinations, namely 

m B Ti = m B F +\f 1 -F 2 - — (F 3 - F A + Ma s {^)F%) , 
^ Vh 

m B T u = F x - 2 ^h(2- u H ) -I) ^ + 2_ _ ^ )F « _ 47ras (/i,)F 6 ") , (15) 
y H {l-u H ) y H 
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where we have suppressed the argument m B u B . Specifically, 

1 

? 



W 1 = X -Tj 



i-u H (i-u H y 

Vh - u H (y H - u H y 

W 3 = . (16) 

2m B (y H - u H ) 
Nevertheless, taking integrals of the form 

1 ni- u H d 3 T u 

dy H dx H K u (x H ,yH,u H )-: — -j—, — , (17) 
u H Ji-yn dx H dy H du H 

with suitable choices of the weight function K u (xH,yH,UH), we can isolate the following 
four linearly independent combinations of the Ff 

(4 - 2u H )m B F + Fi , (18a) 
(1 - u H )m B F + F 2 , (18b) 
F 3 -F 4 + 8na a (tu)Fg , (18c) 

m B F - l -F 3 - X -F± + 47ra s (m)F 6 a . (18d) 

(Recall that we can apply Eq. (14) so that the first combination involves only the leading- 
order shape function.) Here, the treatment of the % dependence in the rate requires care. 
Expanding Eq. (13) in u B ~ A 2 when obtaining the weight function will typically result in 
excessively large coefficients in the u h Fi_ 6 (i71 b uh) terms (which are formally of order A 4 ). 
For example, choosing K u (x B ,yH) = —21x B + 2ly~H + ^5xhVh — ^Vhi we obtain 

1 ff d 3 T u 1 

— // dx H dy H K u (x H ,y H )- — — — — = (1 - 7u H )m B F(m B u H ) + -Fi(m B u H ) + C(A 4 ) , 
T^JJ dx H dy H du H 4 

(19) 

so this eliminates all but the leading-order shape function up to 0(X 4 ) corrections. However, 
we then have the additional contributions 

5 „ , v 49 „ , , 109 „ . , 57 „ . 
-u H Fi{m B u H ) - —u H F 2 {m B UH) ■^u H F 3 (m B u H ) + —u H F 4 {m B u H ) 

83 

- —u H x 4ira s (fii)F^(m B u H ) + 13u H x 4na s (ni)Fg(m B u H ) . (20) 

For this reason, when calculating K u , we keep the full dependence on u B in the rate, rather 
than dropping terms that are formally sub leading in a strict SCET expansion in un/yu ~ A 2 . 
(The analysis of m^-cut effects in B — > X s £ + £" [24, 25] also retained the full u B dependence, 
since doing so facilitates making contact with the total rate in the local OPE [39, 40, 41].) 
Thus, subleading shape functions are eliminated to all orders in u B , and the issue is resolved. 
One straightforward method for obtaining K u (x B , yH, Uh) is then to take different moments 
of the rate with respect to x H and y H , and solve the resulting set of linear equations in the Fj. 
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In Eq. (17), we consider the case where a cut is imposed on p^, i.e. p\ < m 2 D /mB- Different 
or additional cuts will change the limits of integration, calling for different weight functions. 
Table I lists several examples of K u s that isolate the combination itlbF + i*i/(4 — 2uh), 
while Tables II and III give examples that result in (18b) and (18c) respectively. 



(1) 



(3) 



iq 



(2) K? 



1 



1 



9 ( 2 - Uh ) (1 - u H ) 



10(7 -u H )(l-u H )(4 + 3u H )y H 



(454 + 247u H - 71u 2 H )x H y H - 4(1 - u H )(W9 - Au H )y 2 H 



+ 105(7 - u H )x H y 2 H 
5 1 1 



-10(7 - u H )(l - u H )(3A - 27u H )y H 



K? = 



32 (2-u H ) (l-u H ) 8 
+ 2(1 - m^)(2759 - U9u H )x H y H - 525(7 - u H )x 2 H y H 

+ 2{l-u H ){Ul-lZlu H )y 2 H 
15 1 1 



41 {2-u H ){l-u H f 



■2(l-u H y{2W-29u H )y H 



+ (1426 - 1793m h + 157u 2 H )x H y H - 10(109 - Au H )x 2 H y H 



+ (341 - l31u H )x H f H } 



TABLE I: Some choices of K u (xh,Vh,uh) for which the weighted integral Eq. (17) equals msF + 
Fi/(4 - 2u H ). 



(A) Kl 



(B) 



9(i-u H y 



-2(1 - u H )(7 - 15u H )y H - (34 + 7lu H )x H y H 



- 16(1 - u H )y 2 H + l05x H y 2 H 

' 2(1 - u H )(266 - 135u H )y H + 898(1 - u H )x H y H 



(C) Kl = 



Z2{l-u H y 

- 525x 2 H y H + 262(1 - u H )y 2 H 
15 1 



41(1-mh) 7 

- 4:0x 2 H y H + l31x H y 2 H 



58(1 - uufyu + (26 - 157u H )x H y H 



TABLE II: Some choices of K u (xh,Vh,'Uh) for which the weighted integral Eq. (17) equals (1- 
u H )m B F + F 2 . 
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(a) Kl = 



10 



9 1-Uff 



-2(1 - u#)(58 + 32u H + 15u 2 H )y H 



(b) X" 



+ (158 + 104w^ + 53u 2 H )x H y H + (1 - « H )(149 + Qlu H )y 2 H 
- 105(2 + u H )x H f H 

j) ^ r 2(l-^)(92-12^-45^)y H 



16 (l-^) 8 

- 2(1 - Utf)(246 + 139u ff )x H j/H + 175(2 + u H )x 2 H y H 
-2{l-u H ){l& + l7u H )y 2 H 



(c) 



15 



41(1-M 8 



2(1 -u H ) 2 (166 + 93utf)j/ff 



2(483 - 320wf/ - 268u 2 H )x H y H + 5(149 + Qlu H )x 2 H y H 
6(18 + l7u H )x H y 2 H 



TABLE III: Some choices of K u (xh,Vh,uh) for which the weighted integral Eq. (17) equals F3 — 
F 4 + 2F 5 ". 

Now, the subleading shape functions F 5j6 depend upon the light-quark flavour (see Ap- 
pendix A). We indicate this difference between the F 5;6 's appearing in B — > X u tv and 
5 — > X s 7 by using the superscripts 'u' and 's'. In order to cancel the F§ contribution to 
the latter decay, 2 we can use approximate SU (3) flavour symmetry, namely the fact that 



F 5 M - F 5 S m s 



A 



QCD 



(21) 



is suppressed. This enables us to relate the semileptonic process to the radiative process 
and thereby derive an expression for Tq, or equivalently to subleading order. We can 
write 



1 



Kl-Kl + pK^ 



d 3 T u 



dxudynduH 



dx H dy H 



(22) 



m B F(m B u H ) - - [Fi(m B u H ) - 2F 2 {m B u H )] 
- [F 3 (m B u H ) - F 4 (m B u H ) + 2Fg(m B u H )] , 



2 The authors of Refs. [33, 34] have used model-dependent arguments to estimate that the effects of fe.e, 
when integrated over a sufficiently large region, are comparatively small (~ 5%), but that they may cause 
large corrections in the dF/dp^ spectrum for p\ < 0.5 GeV. We avoid any need to consider the reliability 
of these numerics by simply eliminating /s.6, along with the other tree-level shape functions. 
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where 



p(uh) 



(2-^)K~2w) 



(23) 



and F£ = 47ro; s (/ij)F 5 ". If", Iff and if J are any weight functions that give the combinations 
m B F + F 1 /(A—2u H ), (l—u H )m B F + F 2 and F 3 — F 4 + 2F 5 M respectively (examples of which 
are presented in Tables I, II and III). The shape functions in Eq. (22) appear in the same 
linear combination as in the rate dT s /dun- Hence, at NLO we obtain 



1 

ru 
n 



Kl-Kl + pK? 



d 3 T u 



. dxudyuduu 
More generally, we can construct if" such that 



dx H dy H 



1 dY s 



AT 



1 



-AT 



1 

1 



K u (x H ,yH,u H ) 



d 3 T u 



'i-u H y r s du H ' 



-dx H dy H 



(24) 



dxudyuduu 
m B F(m B u H ) + ^(u^F^msUH) 

+ k 2 (u h ) [F 2 {m B u H ) - F 3 (m B u H ) + F A (m B u H ) - 2F^(m B u H )] 



(25) 



For example, we can use 



X u — K u — k u K u 



(26) 



where if^ is a weight function that gives the linear combination m B F + k\F\ + k, 2 F 2 , 
examples of which are given in Table IV in Appendix C (with the corresponding values of 
re" and n 2 shown there). We can also use 



K u = /3K? + 



1-/3 
1-uh 



(iff - if 



(27) 



with P an arbitrary real number (in which case re" = (3/ (A — 2uh) and k 2 = (1— — uh))- 
For any such if", we have 



M" + kI(1-u h )~ 3 M 



+ 0(a s ,A 4 ), 



where M s = (1/T S )M S = (l/T s )(dT s /du H ), i.e. combining M u and M s in this way gives 
an expression dependent only on the leading-order shape function. Taking the ratio of two 
such expressions (two choices of if") at u B ^ then provides us with a relation independent 
of both leading and subleading shape functions. We shall use the superscripts (i) and (ii) 
when we need to distinguish between quantities in the two expressions. We then obtain 



1 

1 o 





"&(«) M «(i) _ &W M «(«)" 






Aii) u(i) Ai) u(ii) 
u ^2 u ™2 


(l-u H )- 3 M s 



(29) 
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where 



6 = (1 - 4) + ( — - u H J (2k? + 4) . (30) 
\m B J 

Since the right-hand side of Eq. (29) is measurable, it enables an experimental determination 
of the CKM ratio on the left-hand side. Additionally, the factor |V^V t *| in this ratio can be 
eliminated by normalizing the photon spectrum by the total B — > X s 7 rate, which is given 
in a local OPE. 

There will be loop and power (A 4 -suppressed) corrections to the rates and hence also to 
Eq. (29). While these are not fully known, one can show that the corrections to Eq. (29) 
are proportional to 

lit) u(i) 
°° , °0 K 2 .... / 31 x 

(multiplied by a s or A 4 ). This needs to be taken into account when selecting {K u< - 1 \ K u ^}\ 
one should avoid pairs of weight functions that result in Eq. (31) being excessively large, 
lest parametrically suppressed terms acquire excessively large coefficients. For example, one 
appropriate choice is to use Eq. (27) for both K u s, with (3^ = 1 and (3^ = 0, after which 
the magnitude of Eq. (31) is less than 1/6 for < uh < m 2 D /m 2 B . 



B. Relations involving B — ► X s £ + £ 

We can also try to isolate shape functions in the process B — > X s l + l~ by taking integrals 
of the form 

iax rl-UH dP"Y U 

dy H I dx H K u (x H ,y H ,u H )- — — — , (32) 

Ji-y H dx H dy H du H 

where 

^max(min) 

27min(max) = 1 — ~T, V " (33) 

(1 - U H ) 

Here, yn = <f/rn 2 B and the low-g 2 region corresponds to 1 GeV 2 < q 2 < 6GeV 2 . However, 
determining K u (xH,yH,un) in the straightforward manner described above proves to be 
problematic in practice. Therefore, we resort to another method, which is based on the 
following observation. Under the transformation xh — > x' H = 2 — uh — IIh — xh, we find that 
f, 1 - H dxn = fi - H dx'rr and 

(1-xh- u h ) <-> (x H + Vh-1), 
(2x H + u H + Vh - 2) <-> -(2x H + u H + yn - 2) . 

This symmetry or antisymmetry can be exploited to obtain K . For example, if K changes 
sign under the transformation, then we can see from the triply differential rate, Eq. (5), that 
integration over xh eliminates the W\ and W 2 terms, whereas the W 3 term remains. Now, 
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Eq. (BIO) shows that F 3 , F 4 and F| occur in W 3 in the same linear combination as in the 
B — > X s 7 rate. 

This still leaves the integration over yjj, and if we choose K u (xh, Vh,uh) = (2x B + % + 
y H - 2) K u (y H ,u H ), where K u (y H ,u H ) satisfies 

fSmax 1 

/ dy H {y H -u H f— (2Re[C 10 „C; a ] + Re[C 10a C 9 *J(1 - f H ))K u (y H ,u H ) = 0, (34) 

then all of the subleading shape functions in Eq. (32) appear in the same combination as 
in the B — > X s 7 rate, which can thus be used to eliminate these functions. Table V in 
Appendix C shows several examples of K u of this form. We observe that z = cos 9 = 
(2xh + uh + Vh ~ 2)/ {yn ~uh), where 9 is the angle between the B and £ + in the center-of- 
mass frame of the . This means that a choice of K u oc {2xh + u b + Vh — 2) is equivalent 
to taking moments of the forward-backward asymmetry, 

d 2 A FB f 1 sign(z) d 3 T 3 r 1 d 3 r 

dz ~^ 1?. 17. z 17~J7. ~JZ ■ ( 35 ) 



dyuduu J_ 1 r dyuduudz 2r J_i dyudundz 

Note also that C 9a is a function of q 2 , and hence of yn (see Appendix B), but in the low-g 2 
region \Cg a \ varies by less than ±1% and we take it to be constant. There is no problem 
taking into account the exact dependence, but integrals over regions of yu must then be 
performed numerically. 

Let T S M S = dT s /du H , and let M u = T^M U and T%M U denote the integrals (17) and 
(32) respectively, with weight functions from Tables I and V. Then we obtain 

r« = l+Jt — (36) 

\ + 2{^-u)k\ M u -k%(1-u h )- 3 M 8 ' 

where re? (k% ) is the coefficient of Fi (F 3 ) in M u (M u ). 

More generally, by the same methods, we can find K u and K u such that 

1 Iff d 3 T u 

- r J" - J^yy K^ H ,y H ,u H)l ^^ dxH iy H (37) 

= m B F(m B u H ) + «i(uh)Fi(771b«h) + K 2 l (u H )F 2 (m B u H ) , 
1 1 dV s 



n ~r s du H 

= -(1 - u H f \^m B F{m B u H ) - ^ [Fi{m B u H ) - 2F 2 (m B u H )] 

- [F 3 (m B u H ) - F 4 (m B u H ) + 2Fg(m B u H )] } , 

M u = ^M u = If K u {x H ,y H ,u H ) d T dx H dyu 
r o ™JJ dx H dy H du H 



m B F(m B u H ) + ^k 2 {u h ) [Fi(m B u H ) - 2F 2 {m B u H )] 
+ K,f(u H ) [F 3 (m B u H ) - F A (m B u H ) + 2F£(m B u H )] , 
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where F 5 S = 47ra s (/ij)F 5 s . Tables IV and VI show (further) examples of such weight functions, 
along with the corresponding values of the coefficients re" 2 and k% 3 . Then 



[4 - 4 e ]M u + k u 2 M u - (1 - u H )-*M s 



(38) 



= ( 



;K 



+ 0(a s ,A 4 ), 



so in this case we have a combination of M u , M s and M w that is dependent only on 
the leading-order shape function. Taking the ratio of two such expressions (two choices of 
{K U ,K U }, denoted by superscripts (i) and (ii) as previously) at uh ^ then provides us 
with another relation independent of both leading and subleading shape functions. 
Specifically, let 



M s = —M s = 



if k% = 



If 



[K% - K l i\M u , if K» ^ 

f ««M", if «f ^ K | 



(39) 



if «f 



1 



hu 2 ft 3 



(1 - u H )- 3 M s , if «f ^ k| 



-/4 £ (1 - m^)- 3 M s , if 



AC 



It 



and 



c 



(«« _ K « + K « + + - UH ) («« - «f ) (2«? + «!?), if ± and «« ^ k| 



1 + 2 



if = 

if k 2 = tdf 



We find that 



1 

1 



cfO^uW _ c f M u{u) 



c^\M s ^ + rM*W) - $>(M a M + rM u W) 



CO, 



(40) 



(41) 



where r = I^/T^, or 



1 



W, 



(42) 



In the special case where = and k 



i) 



«" w , Eq. (42) reduces to Eq. (36). 



The loop and power (A 4 -suppressed) corrections to Eq. (41) can be shown to be propor- 
tional to 



eg 



cf [«^(1 + ««)](« 



+ 



(ii)r 



;(0r 



(43) 



+ 
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where Co = (idf — + k% + K^nf) + (A/me — uh)(k^ — i^){2k^ + K V)- When selecting 
{K u ^\ K u ^}, {K u ( ll \ K u ^}, one should avoid those sets of weight functions that result in 
Eq. (43) being excessively large. The following combinations of weight functions are suitable 
choices: 

K u(i) = (^ ( 2 ) or (3) [Table I] and K U W = (7), (8) or (9) [Table V]; 
K u(€) = ( 5 ) Qr ( 6 ) [ Table IV ] ; K u(i) = ( 10 ) ) (ii) or (12) [Table VI], 
and K U W = (7), (8) or (9) [Table V]. 

C. Perturbative Corrections 

Let us now consider the feasibility of incorporating perturbative corrections in our rela- 
tions. In Ref. [32], the complete set of subleading corrections (to all orders in a s ) for the 
triply differential spectrum of B — > X u iv was derived. It was shown that prohibitively many 
new shape functions appear at order a s AQ CD /m b , and hence it is not phenomenologically 
viable to work to that order. 3 However, one may choose to work to order (a s \°, a° s \ 2 ), by 
including perturbative corrections to just the leading-power terms. Recall that there are 
two perturbative scales, /if, ~ mi, (hard) and /ij ~ a/^vV^cd (jet). It is straightforward to 
take into account the relevant hard corrections. Including the effect of corrections to the jet 
function which is convoluted with the shape function /(°), is more involved: one has to 
"invert" a distribution (see Eq. (8)). An implementation akin to Refs. [26, 27, 28, 29] is left 
for future work. Nevertheless, before this is done, we can still use the less direct approach 
mentioned in the introduction, using two instances of Eq. (28) or (38), with appropriately 
modified right-hand sides. For example, one can extract the leading-order shape function 
from the analogue of Eq. (38), with K from Table V, and substitute this function into a 
second choice, with K u from Table I. Finally, we note that the extent to which Eq. (29) or 
(42) varies with respect to u H or different combinations of the K u s and K u s will provide a 
measure of the effect of a s and A 4 corrections. 

IV. CONCLUSION 

In this paper, we have established a method for obtaining |V^b]/|Vt&V^*| that includes 
O '(Aqcd / mb) corrections in a model-independent way. Our approach relies upon a class of 
relations between the inclusive decays B — > X u tv and B — > X s ^ that are valid including 
the first-order power corrections (see Eqs. (24) and (29)). Alternatively, one can use a 
separate class of relations involving B — > X s £ + £~ (see Eqs. (36) and (42)). Experimentally 
required cuts make shape-function effects important in these processes. Their differential 
decay spectra in the shape function region have previously been derived to subleading order 
with the help of the Soft-Collinear Effective Theory. These rates involve a number of non- 
perturbative but universal shape functions in different linear combinations. We are able 
to eliminate these sources of hadronic uncertainty by taking suitable weighted integrals of 
the triply differential rates. Hence, our results incorporate NLO power corrections while 



3 Unless these shape functions appear in the rates in only a much smaller number of linear combinations. 
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avoiding model dependence. There are many possible weight functions (see e.g. Eqs. (26) 
and (27)); different choices provide a consistency check on the determination of \V u b\. 
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APPENDIX A: SHAPE FUNCTIONS 

The leading-order shape function is 

= l - (B v \h v 5(£ + -tn-D)h v \B v ) , 
where h v is the heavy quark field. The subleading shape functions are 



(B v \O (£ + ) 


B v ) 


= /o (2) (^) 5 


(Bv\0^t + ) 


B v ) 


V n-v / 


(B V \P?\£ + ) 


B v ) 




(B v \Of(£t 2 ) 


B v ) 






B v ) 


V n-v/ 


n a np(B v \Of(£f^ 3 ) 


B v ) 




i^){B v \Of{£t,2,z) 


B v ) 





(Al) 
(A2) 



where = g» v - (l/2)(rj¥ + n v n») and = {l/2)e^ uaP n a n (3 . The ultrasoft operators 
are 

r dx~ 



of (it, 4) 



8tt 



_i T -/+ 

e 2 



(A3) 



d 4 yT [h v (x)Y(x,0)h v (0) iO h (y)] , 
h v {iDi,S(e + -in-D us )}h v , 
h v [iD^,5(£ + -in-D us )]^ l5 h v , 
h v 5(£t-in-D us ) {iDi^iD^} 5(£+ -in-D us )h v , 
-- h v b[£\ -in- D us )gG^ ± 5(£l[ -in- D us )^ b h v , 

\ {h v 5{£t-in-D us )^P L T A q n } 5(4 - in-d) {f^P^t -in-D us )T A h v ) , 



where x^ = n-xn M /2. Here, Oh is the NLO term in the HQET Lagrangian, Y is an ultrasoft 
Wilson line, igG^ sA _ = [iD^,iD^] and q™ s = ($if)/4q us . The operator Of, which appears 
in the definitions of / 5j6 , depends upon the light-quark flavour, u or s. 
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APPENDIX B: HARD COEFFICIENTS 



In this appendix, we present expressions for the hard coefficients in B — > X u £u and 
B -»• X s £ + ^ [24, 32, 36]. At lowest order, we have 

rPx 

wf> = hifa,m b ,ii) / dk + j(°\p- x k + ,fi)f(°\k + +A-p+,ri. (Bl) 

Jo 

For _B — > X u £u, we have 

fc? = ^[Cf>] 2 , (B2) 

{Vh-uh) 4 (y H -u H ) 2 

3 2m B (y H -u H ) ' 



where 



Cf>(*,l) = l-^^{21n^) + 2Li 2 (l-^) + ln(^)(^) + ^ + 6 
1) = 



1N _ a s {m b )C F \ 2 { 2u>ln(a;) 



4tt \(1 -w) (1 — cD)2 J ' 



1N _ a s (m b )C F \ (1 - 2u>)u> ln(c2>) cD 



Here, a) = uo/m b . 

For 5 -> X s £+£-, we have 

/,« Vlrl 2 ^lr |2\ 2Re[C 7 C 9 *] 2|C 7 | 2 

2(l—u H ){\ n a \ n i2 r - i\ Ciob 8 C 7 



h 



(Vh-uh)^ V 2 (1-Vh)(vh-uh) 

-4Re[C 10a C 7 *] 2Re[C 10a C 9 *] 



''3 



m B (l - y H )(y H - u H ) m B (y H - u H ) 



The full expressions for the coefficients C 7 ,9,ioa,iof> are given in Ref. [24]. When we ignore 
0(a s (m b )) corrections, they simplify to 

(B5) 

c 7 NDR M , 



c 9 


— C 9 a — 


/nmix 
°9 ' 


c 7 


= c 7a = 


m b {no) 






m B 


ClOa 


= ClOa = 


Cio ) 


Ciob 


= o, 
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where /i ~ mj and 



C 9 mix (/i ) = C 9 NDR (/i ) + ^ (3C 3 + C 4 + 3C 5 + C 6 ) - ^(1, s) (4C 3 + 4C 4 + 3C 5 + C 6 ) 

777 I 

+ , s) (3d + C 2 + 3C 3 + C 4 + 3C 5 + C 6 ) - -/i(0, s) (C 3 + 3C 4 ) 
+ 0(a fl (A*o)) • (B6) 
The function s) is given by 

= iM£>-|m* + £ + !c-§<2 + 0>/iI=Ci 



X 



0(1 - C) ( - m + In 1 + ^ | ) + 9(C - 1) 2 arctan 



h ^ = h + l^-^ 8+ ^ (B7) 

with £ = 4z 2 /s and s = g 2 jm\. 

In the expressions above, Ci_6, C^g 311 , Cio are the coefficients of the corresponding op- 
erators in the effective Hamiltonian for b — > s£ + £~ (for which the NLL calculations were 
done in Refs. [42, 43]), while C«^ nix differs from Cf of Ref. [42] by only an 0(a s ) piece. 
Note that there is a complication in the perturbative power counting. Above the scale nib, 
one usually expands in a s , with a s \og(mw/rrib) = 0(1). Because of mixing with Oi t 2, 
C 9 ~ log(m^/m&) ~ l/tt s , whereas C 7jW ~ 1. However, numerically |Cg(m&)| ~ Ci . 
This problem is exacerbated by the fact that in the shape function region only the rate is 
calculable, not the amplitude. The solution is to use a "split matching" procedure, which 
decouples the scale dependence above and below /i = and thereby allows us to consider 
the coefficients as 0{1) numbers in the latter region [24]. 

At next-to-leading order, we have 

Zrrib J 

rPx 

/ dk + J^\n-pk + ,fi) fj: 2) (k + + r + ,fi) 
Jo 

+ ]T^^ fdkidktJ i *- 2 \n- V k+^)ff\k+ + r\^) 

+ fdktdktdkt J^\n-pk+,rif<?\k+ + r + ,VL) 

n-p J 

+ ..., (B8) 

where r + = A — p^-, j = 1,2 and f = 1,2,3. The ellipses denote terms that have jet 
functions J that start at one-loop order or higher. (These terms are given in Ref. [32].) 



+ V- I >'i f ("■!>) 



4 



r=l ' "° 







K f (n 


p) 


m b 




h r /(n 


p) 
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When we keep the full dependence on uh, the /i™ 3 are 
h m_ 1 h m_ (l-u H )(2-y H -u H ) 

n l ~ o ' 2 ~~ o/- \9 ' 1 



I ii _]_ h m = (I - u H )(2 - y H - u H ) hlu = 1 

8 ' 2 2(y H - u H ) 2 ' ^m B (y H - u H ) ' 



2„_ 1 U 2u _ (l-u H )((4-u H )y H -f H -2) 2u _ 1 



1 4' 2 yH{ljH-u H ) 2 ' 3 2m B (y H -u H y 

h 3u = l _ h 3u = ( l ~ U h) h 3u = 1 

1 ^Vh ' 2 yH{yH-u H )' "' 2m B y H (y H - u H ) : 



% Au 1 = (1 -Mj/)(2-t/ g -M g ) ^ 4u = 1 

1 4f/# ' 2 Vh^Vh-uh) 1 2m B y H (y H - u H ) : 



,5«__i ,sn _ 2(1 -n H )(l 5n _ 1 

1 — o ' ""2 ~ I- „. \9. ' 3 



2' " {yH-u H ) 2 m B (y H -u H )' 

2(1 -u H ) 2 
(Vh ~ u H ) 2 ' 



h^ = 0, = , 2 _ (1 Uh) L , h 6 3 u = 0, 



and the /i^ 3 are 

4|C 7a | 2 -(|C 10a | 2 + |C 9a | 2 )(l-^ 



4(1 -y H ) 2 

(2 - j/ H - (4|C 7a | 2 - (|C 10a | 2 + |C 9a | 2 )(l - y H ){\ - uh)) 



h 3 



(1 -Vh){vh -u h ) 2 

id Re[C Wa Cg a ] 



2U 



K 

h 2U 



m B (y H - u H ) ' 
4|C 7a | 2 -(|C 10a | 2 + |C 9a | 2 )(l-^ 
2(1 -Vh) 2 



4|C 7a | 2 2 \ H VHUH + 4Re[C 7a C 9 *J(2 - jfo - Ujr ) 
1 — 



- u H ) 

+ (|Cioa| 2 + |C 9 „| 2 )(2 - 4y H + f H + Vhu h )(1 - u H ) 



y^tt _ 2Re[Ci0a Cg* a 



/if 



/if = 2 



TUB (l/H ~ U H ) ' 

4|C 7a | 2 + 4Re[C 7a C 9 * a ](l - gg) + (|C 10a | 2 + |C 9a | 2 )(l - y H f 

2////H - !J„) 2 
MC 7a \ 2 - (\C 10a \ 2 + |C 9a | 2 )(l - y H )(l - u H ) 



Vh(1 - u h )(uh - u H ) 



h , u = 2 2Re[C 10a C* a ] + Re[C 10a C 9 *J (1 - y H ) 



m B y H {l - yn)(yH - u H ) 
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hf = 



h 4U 
u 3 



h\ li 



hf 

hf 

hf 



A\C 7a \ 2 + ARe[C 7a C 9 *J(1 - y H ) + (|C 10a | 2 + |C 9a | 2 )(l - y H f 



ynijjH - u H ) 2 



4|G 



2y H (l -y H ) 2 
, 2 2-y H -u H 



7a | 



+ 8Re[C 7a C* a ](l-u H ) 



= -2 



= 4 



i - yH 

+ (\C 10a \ 2 + \C 9a \ 2 )(2 -y H - u H )(l - u H ) 

2Re[C 10Q C 7 * a ] + Re[C 10a QJ (1 - gg) 
■m B y H {\ - y H ){y~H - u H ) 

A\C 7a \ 2 + 4Re[C 7a QJ(1 - y H ) + (|C 10a | 2 + |C 9a | 2 )(l - y H f 

(i-m) 2 

(4|C 7a | 2 + 4Re[C 7a C 9 *J(1 - y H ) + (|C 10a | 2 + |C 9a | 2 )(l - y H f) (1 - t/g) 

(l - y~H){y~H - u H ) 2 

2Re[C Wa C 7 *J + Re[C 10a C 9 *J(1 - y H ) 



m B (l- yn)(yn ~ u H ) 



= 0, 



= 4 



A\C 7a \ 2 + 4Re[C 7a C 9 *J(1 - u H ) + (|C 10a | 2 + |C 9a | 2 )(l - u H f 



(yH - u H ) 2 



0. 



APPENDIX C: WEIGHT FUNCTIONS 



(4) 



(5) 



(6) 



-"iv 



1 (2x^ + ^ + ^-2) (1 -y H ) (2yH-u H - 1) 



(yH - u H ) 



(1 + uh- y^ 



N(u H ) 

N M = 4(! - «ff) 2 (l " 14u h " 94«h - Uu% + ujj) - 2*4(1 - O logMH 



30 
1 

2 ' 



k 2 



-1 



«i = 



6 [ - (7 - + 4x H !/ g + 6y^] 

7 -u H 



10' '" z 5(1 -u H ) 
105 [14(1 - un) 2 yH + 5(2 + 7u H )x H y H - A5x H y H ] 



101 

2 

101 ' 



(i-uff) 6 

109 -4mh 
101(1 



TABLE IV: Some choices of K u (xh, Vh,uh) for which the weighted integral Eq. (17) depends only 
on the shape functions F, F\ and F 2 . The coefficients k\(uh) and k 2 (uh) are defined in Eq. (37). 



19 



/ ? x R u = 1 (2x H + u H + y H -2) (y m -y H -u H ) 2 
N(u H ) (y H - u H ) (Vh* ~ Vh) 

x{A- B[l - (y H *-y H ) 2 }} (Vh* - 2y H ) 

N(u H ) = 8(1 - u H ) [ dy H (y H -u H ) 2 ^ H * ^ H u ^ (y m _ 2y H ) 
Jy min [Vh* - Vh) 

x{A- B(l - y H )} {A - B[l - (y H * -y H ) 2 } } 
u 8(1 -u H ) ft™ _ Ay H *-y H -u H ) 



(9) K 



N(u H ) 



r max w , - n {yH*-y H -u H V _ v 

/ dy H {yH-u H ) (1 - yjf) — 73 r-^ — (y#* - 2y#) 

J Smin tJ/H* - yn) 

xB{A-B[l-(y m -y H ) 2 ]} 



(8) = N ( Uff j ( 2x h + u H + y H - 2)( y y H *-y H -u H ) 3 y H ( y yH* - 2y H ) 



x {A-B[l-(y H .-y H ) 2 ]} 



"J/max 

AT(u^) = 8(1 - u H ) I dy H (yH -u H ) 3 (yH*-y~H -u H ) 3 y H (y m - 2y H ) 



j. i .i j. 1. 1 



x{A- B[l - (y H *-y H ) 2 }} {A - B{\ - y H )} 

) [ l 

N(u H ) 



K 2 = ^77 — v~ I dy H (yH-u H ) 3 (l - yH){ym-yH-u H fy H (y H * - 2y H ) 



mm 



xB{A-B[l-{y m -y H ) 2 ]} 
1 (2x H + u H + Vh - 2) (yH*-yH-u H ) 2 (Vh* - 2y H ) 

' N{u H ) (y H -u H ) (y H *-y H ) {A-B(i-f H )} 

N(u H ) = 8(1 - u H ) [ dy H (y H -u H ) 2 ^^—^— -^—{y H * - 2y H ) 
Jy min (Vh* ~ Vh) 

{A-B(l-y H )} 
{A-B(1-? H )} 
u 8(1 -u H ) _ (y H ,-y H -u H ) 2 

^ = —R^L *«<»»-"'> f 1 - »») {Sh ,- Sh) (»- - 

X 



M-B(l-^r)} 



TABLE V: Some choices of K u (xh,Uh,uh) for which k^(«h) = k<i(uh) m Eq. (37). a Here, 
.4 = -2Re[Cioa C 7 *J, 6 = Re[Ci 0a C 9 * a ] and = £/ min + y max . 

a Note that Example (9) requires a harsher cut, e.g. 2GeV 2 < q 2 < 6GeV 2 (rather than 1 GcV 2 < q 2 < 
6GeV 2 ), so that it is not singular. 
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(10) 




1 (2x H + UH + VH-2) 
N(u H ) 4(1 -u H ) 




N{u„) 


= 2 / dy H {y H -u H f {A - B(l - y H )} 




r"2 










(11) 


K u 


1 (2x ff + u H + y H - 2)y H 
N{u H ) 4(1 - UjT ) 




N(u H ) 


= 2 / dj/ ff j/H(j/ff-«H) 3 M-^(l-j/H)} 

"'j/min 




hu 2 


= N , \ / <*j/tfB(l ^) 3 y^ 

^ ' 2/min 




K U 
o 


2 rVm&x 

= ; — r / dyu(yu— Uff) 3 {A — B(l — Vh)} 


(12) 


K u 


1 (2x H + u H + Vh - 2) (yH*-yH-u H ) 2 _ s 
7V(?j h ) 4(1 - u H )(yH - u H ) (Vh* - Vh) h * 




N{u H ) 


Jy min (Vh* ~ Vh) 

x{A- B{1 - y H )} {A - B(l -y H *+y H )} 




hu 2 


i^j Jy miri {yH* - Vh) 
x (y H ,-2y H )B{A-B(l-yH*+y H )} 






= 



TABLE VI: Some choices of if and the corresponding coefficients k*[{uh) an d k z{ u h), which 
are defined in Eq. (37). Here, A = -2Re[Ci 0a CfJ, 6 = Re[Ci 0a C 9 * a ] and = y min + y mgx . C 9a 
may be taken to be constant, in which case the integrals can be evaluated analytically. 
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